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ABSTRACT 

Fundamental solution of Fokker - Planck equation is built by means of the Fourier 
transform method. The result is checked by direct calculation. 
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We see from recent publications (ref. [1-5]) that the Fokker - Planck equation is an object of steady 
investigations. It is a very useful modelling tool in many fields. Therefore it is desirable to have exact 
expression for fundamental solution. 

The object of our considerations is a special case of Fokker - Planck equation, which describes evolu- 
tion of 3D continuum of non-interacting particles imbedded in a dense medium without outer forces. The 
interaction between particles and medium causes combined diffusion in physical space and velocities space. 
The only force, which acts on particles, is damping force proportional to velocity. 

For this special case the Fokker - Planck equation is: 

dn dn dn dn f d d d \ fd^n d^n d^n\ 

di^'^rx^'^d^^'^Tz-''[d^,^'^''^^d^/'^''^^d7^^^^ 

where 

n — n{t, X, y, z, Vx^y^z) ' density; 

t - time variable; 

x,y, z - space coordinates; 

v^VyV^ - velocities; 

or - coefficient of damping; 

k - coefficient of diffiision. 

We shall search solution of equation (1) for unlimited space 

— oo < X < +00, — oo < y < +00, — OO < Z < +00, (2) 
—oo < < +00, — oo < Vy < +00, — oo < < +00. 



Let us denote by mq initial density 
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noix, y, z, v^VyV^) = n(0, x, y, z, v^VyV,). (3) 
We shall use the Fourier transform method. Let us denote by A'^ Fourier transform of density 

N = Nit, Pjc,Py,Pz,1x.(ly,1z) = (4) 



OO OO OO OO OO OO 



Iff \ \ \ \ \ ^ expi-iixp^ + ypy + zPz + V;,q;, + Vyqy + v^q,))ndxdydzdv^dVydv,. 



—OO —OO —OO —OO —OO —OO 



and by Nq - Fourier transform of initial density: 

NoiPx,Py,Pz,9x,9y,9z) ^ "^(0. Px,Py,Pz,^x,9y,9z) (5) 

where Px,Py,Pz " space coordinates momentum variables; 

^x,^y,1z ' velocities momentum variables. 

Multiplying (1) by e.x^{-i{xpx + ypy + zp^ + v^q^ + Vyqy + v^qj) and integrating over whole space 
and velocities, we obtain 



dN dN dN 3Af I 2 2 2 , 

+ (ocq^ ~ 9^ ^"'^y ~Py'dq~^ ^"'^^ ~^^^dq~^ y'" ^"/^' 

The equation (6) is a linear differential equation of first order, so it can be solved by the method of 
characteristics. Relations on the characteristics are 

dpx dpy dp^ dq^ dqy dq^ dNIN ^ 

aq^-Px aqy - Py aq^ - p, k {ql + ql + ql) ' 

First three equations (7) have three integrals: 

Px - const; (8) 
Py = const; 
p^ = const. 

If we combine (8) with next three equations (7), we get three further integrals 

(aqx - Px) c~"' = const; (9) 

(aqy - Py) e~"' = const; 

{aq^ — p^) e~"' — const. 



We solve (9) for qx, qy, q^ and obtain 



Px , at { Px I /I /A\ 

qx^— + e"'\ qxo ; (10) 



Py . at Py 
qy^— + e"'\ qyo - — 
a \ a 



n - A- in - 
1z- + ^ \1zO 

a \ a 
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On the other hand, if we solve (9) for q^Q, qyQ, q^Q, we obtain 



Px ^ -at Px 

a \ a 



(11) 



Py -a, Py 
a \ a 



^zo - — + e \q^ 

a \ a 



To get the last integral of (7), we replace current velocities momentum variables q^, qy, q^ by initial veloci- 
ties momentum variables q^^y, qyo, q^o inq^ + q^ + q^ 



ql + q] + ql = -^\pl + p] + pl | + 



(12) 



+ ^ £"'[ p^iq^o - ^) + PyilyO - ^) + PziizO " ^) 1 + 



a 



a 



+e""\(q^o-^f + (qyO-^f + (qzO-^f 
^ a a a 



Integrating (12) in t, we obtain the last integral 

IniN) + k 



{pl + pl + pl) + 



(13) 



+ ^"'\ P^^^^o - ^) + PyilyO - ^) + PziizO - ^) 1 + 



Py 

a 



Pa 
a 



y"'i< p^ 

+ (1x0 

2a V oc 



? + ilyO-^f + ilzO-^f 

a a 



const. 



We replace in (13) initial values of velocities momentum variables by their current values 



ln{N) + k 



(pI + pI + pI) + 



(14) 



+ I Pxilx - — ) + Pydy - — ) + Pzdz - — ) 1 + 

I a a a 



1 

2a 



+ ^\ilx- + ily - + ilz - ^) 



a 



Pz-,2 

a 



= const. 



To determine the constant term in (14), we write the same expression for initial values and equate both 

expressions 

ln(N) + k\^(pl + pl + pl)+ (15) 



2 ( Py ^ 

+ Pxilx ) + Pyily ) + Pzilz ) + 

\ a a a ) 
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+ — \iqx-—) +i9y-^) ) 



2a 



a 



a 



a 



IniNo) + k 



a 



a 



+ 2^ ( (9.0 - + (9,0 - + (9^0 - 



a 



Solve (15) for ATq 



a a 



■ exp -k 



4 (^' + Z'' + P') + A (1 - e""0L(9. - — ) + ^,(9, - — ) + pMz - -) 1+ (16) 



+ J_ (1 _ eM{q. - + (9. - -f + (9. - -f 
2a \ a a a 



The A^o(- ■ ■ ) in the (16) means, that one have to calculate from initial density according to (5) and 
then replace values of its arguments by expressions (11). 

Let us specify initial density value as product of delta functions 

fiQ = S(x - XQ)S(y - yo)3{z - Zq)S{v^ - XQ)S{Vy - Xq)S{v^ - xq). (17) 

The Fourier transform of initial density (17) is 

No ^ expi-iixop^ + yoPy + zoP^ + v^^q^ + v^o^y + v^^qj). (18) 



Substituting (11) for Nq arguments in (18) gives 

1 \ 



exp < -I 



XoPx + yoPy + ZoPz + 



(19) 



a 



+v4 ^ + e-"'iq, - ^) 1+ v,o( ^ + e-"\qy - ^) )+ v,of ^ + e-(^, - ^) 



a 



or 



1 

(2T)6 



p,\ xo + ^ (1 - e-«') + Pyho + -f (1 - ^ ) + Pz\zo + ^ (1 - e ) • (20) 



+e "'{qxVxO + qyVyo + q^Vzo) 
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It is clear that (20) is the Fourier transform of (21) 

no = s(x - (xo + ^ (1 - e-"'))] sfy - (yo + ^ (1 - g-))] sfz " (^o + ^ d " e""'))] (21) 



51 V, - v^oe-"' \S\Vy- Vy^e-"' \d\v,- v.ne""' I- 



To get inverse Fourier transform of n from its Fourier transform (16) we use two known results (ref. 



[6]): 



1. A product of 2 functions A(co)B(co) tranforms to convolution a(x) * b(x) (where n is a number of 

(2;r)" 



independent variables). 

2. The Gaussian exponent of quadratic form e~"'' with matrix A transforms to exponent of quadratic 

\ det(/ 
In our case the matrix is 



form with inverse matrix a / e 4 

det(A) 



A = k 



2a^ 



a'- 



1 

2'^2 



(1 - e-"') -^{l- e-^"') 



la 



(1 - e-^"') 



(22) 



The determinant is equal to 



^at{\-e-^'")-2{\-e-"') 



al\l 



det{A) = k 



(23) 



Let us denote by D expression 



-at\L 



D = 



det{A) _at{\- e'^"') - 2 (1 - e'"') 
2^4 



(24) 



The inverse matrix is 
1 

W 



2a 



-2al 



2a^ 

i (1 - e-') -id- e-^-) ^ - ^ (1 - .-) . i (1 - .-) 



(25) 



Combining (21) with (25) we obtain expression for the fundamental solution: 

n3 



1 

2a 



(l-e-2«') x' + / + z' 



(26) 



■ (1 - - ^ d - e"'"') 1 1 XV, + yv, + zv, 1 + 



a- 
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where * means convolution of two functions. 



The convolution of arbittary function with product of delta functions simplifies to substitution delta 
function arguments for this function arguments. Finally, we obtain 



1 

la 



(27) 



(1 - e-«0 - ^ (1 - I Uv, + yv, + zv J + 



where 

X = X - (xo + — (1 - e-"')); = - (>'o + — (1 - z = z - (zo + — (1 - e-"')); (28) 

a a a 

Vx^v^- v^oe'"'; Vx = Vy- Vyoe'"'; v^^v^- v^oe""'. 

This is the fundamental solution of Fokker - Planck equation. 

Let us check vahdity of solution (27-28). Direct differentiation of (27-28) and substitution to (1) 
leads to cumbersome calculations. Therefore we use "senu-reverse" method. (27) has Gaussian form, so 
we search Gaussian solutions of (1): 



n = exp ] A(f)\ x^ + y^ + \+2B{t)\ xv, + yv, + zv J + C{t)\ vl + v] + v\\ + 3Qit) \. (29) 



To get rid of exponents we write 

/ - /n(n) - Ait{x^ + y^ + z^ |+ 2B(0| xv^ + yVy + zv,\ + C(0[ vl + vl + vl\ + 3Q{t). (30) 



/ must satisfy equation 



3/ 3/ 3/ 3/ 



+ Vr^ — ^v^^ — \-v,^ a\Vr 1- Vv 1- V 



dt 'dx 'dy 'dz 



3/ 3/ 3/ 



dv^ ' ^ dvy ^ 3vj 



■ 3a 



(31) 



3v? 3v?, 3v? 



3/ 



3v, I \dvr I 1 3v 



3/ Y ( 3/ 
+ 



instead of equation (1) for n. 

Substituting (29) for A, B, C, Q in (30) and collecting of similar terms leads to following equtions 

dA 



dt 



■- 4k B^; 



dB 
dt 



= aB-A + Ak B C; 



(32) 
(33) 
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^^2aC-2B + 4kC^; (34) 
at 

do 

-i^^a + lkC. (35) 
at 

It is not easy to solve nonlinear system (32), but our task is simpler. We have only to check, that 

at (\-e-^"')-2(\-e-"'f 



Mt) = (-r^](^il-e-"'')l (37) 



\4kD ) \2a 



D 

1 



^«H^^Jl^-^^'-^""'^"i^^-^""^|^ '''' 

Q{t)=-]^ln{D{t)). (40) 

satisfies equations (32-35). Direct substitution proves this. 
We proved vaUdity of (27) for the special case 

Xq - 0; yo - 0; zq - 0; v^q - 0; v^q - 0; - 0. (41) 



For the common case we prove that differential operators 

a a a 
ax ' ay ' az ' 



(42) 



and 



-(l-e-«0^ + e-"^^; - (1 - e""') ^ + e""' ^ ; - (1 - g-'^O ^ + e""' ^ ; (43) 
or dx av^ a ay aVy a az ov^ 

are symmetries of PDE (1). This statement is trivial for (42) because (1) does not contain x explicitly. To 
prove the statement for (43) we build prolongation of differential operator according to Lie prolongation 
formula (ref. [7]) 

/^1=A(^«")-Sa(^A (44) 



^ax'j ax* 

where 

u" - a set of dependent variables; 
x' - a set of independent variables; 
Di - fuU derivation on x' operator; 

Su", Sx^ - actions of infinitesimal synometry operator on variables. We use this non-standard notation 
instead of usual f , to empasize their nature as small variations of variables. 

6\ -zr— - induced action of infinitesimal symmetry operator on derivatives. 

I ax' j 
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For operators (43) Lie formula gives 

-il-e-'") — -e-"'— + e-'"\-^ + a^\^; (45) 
a ax dv^ \ ax av^Janf 

- (1 - e-"') — - e-"' — + e-"'\ " ^ + « ^ ^ ; (46) 
a ay avy \ ay avyjan, 

a az aVj dz av^jan, 

dn 

where = — . For second order derivatives we have <y(n„„) = 0, ^(nvv) - 0, •^(Iwh') - 0. 
a? 

It is easy to calculate, that the action of (45-47) on (1) is identical zero. We proved, that differential 
operators (42-43) are symmetries of PDE (1). 

The action of operators (42-43) on solutions (41) increases variables XqJq ■ ■ ■ from zero to arbitrary 
values. In this way we get from (41) solutions (27-28). 

We checked the solution (27-28). 



DISCUSSION 

The main result of these considerations is the closed form expression for fundamental solution of Fokker - 
Planck equation without forces. This result can be helpful for modelhng of more sofisticated problems. It 
demonstartes in simple and clear way the result of combined diffiision in physical space and velocities 
space and inertial motion. 
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